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Abstract. We consider saddle point integrals in d variables whose phase func- 
tion is neither real nor purely imaginary. Results analogous to those for Laplace 
(real phase) and Fourier (imaginary phase) integrals hold whenever the phase 
function is analytic and nondegenerate. These results generalize what is well 
known for integrals of Laplace and Fourier type. The method is via contour 
shifting in complex d-space. This work is motivated by applications to asymp- 
totic enumeration. 



1. Introduction 



Integrals of the form 
(1.1) X(A) J(A; (/., A) I e'^'t'^''^ A{x) dx 



arise in many areas of mathematics. There are many variations. This integral can 
arise in one or more variables; the variables may be real or complex; the integral 
may be global or taken in over a small neighborhood or oddly shaped set; varying 
degrees of smoothness may be assumed; varying degrees of degeneracy may be 
allowed near the critical points of the phase function, (f). Often what is sought is a 
leading order estimate of X(A) as the positive real parameter A tends to cxd, or an 
asymptotic series ~ (^n9n{^) where {gn} is a given sequence of elementary 
functions and the expansion is possibly nowhere convergent but satisfies 

N-l 

(1.2) 1(A) - c„g„(A) = 0(gjv(A)) 

n=0 

for any > 1. 

In recent work on the asymptotics of multivariate generating functions [PW02; 
PW04; PW08; BBBP08; BP08; BP04], we have required results of this type in 
the case where the phase function (p and the amplitude function A are analytic 
functions of several variables. The phase function is typically neither real nor 
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purely imaginary, but satisfies Re{0} > 0. Although these results are known, at 
least at a folklore level, we could not find this case analyzed in the literature in 
the generality we required. We therefore wrote this note to patch the gap and to 
remedy what would have been ghost citations in [PW04, Lemmas 4.7 and 4.8]. 

To better explain the relation between the present results and the literature, we 
discuss two special cases in which the results are well known, but by substantially 
different methods. 

Integrals of Fourier type. Let /, A : R'^ ^ R be smooth (that is, C°°) 
functions of d real variables, with A having compact support. Taking </> = —if 
gives the Fourier-type integral 

J(A) = j e'^^^'-'^A{x)dx. 

The standard method of studying this integral is as follows. If / has no critical 
points in the support of A, then integration by parts shows 1(A) to be rapidly 
decreasing: J(A) = 0{\~^) for all positive integers, N. Using a partition of 
unity, the integral may therefore be localized to neighborhoods of the critical points 
of /. At an isolated critical point V/ vanishes; if the Hessian matrix is non- 
degenerate then the Morse lemma produces a smooth change of variables under 
which f{x) = S±{x) := X]j=i ^■'^j- This reduces the gc^neral problem to the case 
(j) = iS. To solve this, expand A in a series. Each term may be explicitly integrated, 
resulting in an expansion in decreasing powers of A: 

(1.3) J(A) = ^c„A-''/2-". 

n>0 

To see that the resulting series for X(A) satisfies (1.2), one uses integration by parts 

again to bound the remainder term. The coefficients {cn} are determined by the 
derivatives of A and (j). A particularly lucid treatment of this may be found in the 
first two sections of [Ste93, Chapter VIII]. 

Integrals of Laplace type. The other well studied case is the Laplace-type 
integral, where (j) is real. Localization of the integral to the minima of (f> is immediate 
because the integrand is exponentially small elsewhere. Integrating over balls whose 
radius has order A~^/^ shows that near a quadratically nondegenerate minimum, 
xq, the bound A{x) = 0(|x — a;o|^) translates into the bound 

(1.4) X(A) = ©(A-^'^+^^Z^) _ 

Again, one changes variables, expands A into a power series, and integrates term 
by term, to obtain the series 

(1.5) J(A) = ^c„A-(''+")/2. 

n>0 

Applying (1.4) to the integral of the remainder term shows that (1.5) is an asymp- 
totic expansion for 2"(A). Classical treatments of integrals of Laplace type may be 
found in many places [BH86; Won89], often accompanied by separate treatments 
of the Fourier case. 
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Complex methods. The series (1.3) and (1.5) are formally identical. Fur- 
ther inspection of formulae for c„ in the literature shows these to be nearly the 
same as well, differing only by constant factors of unit modulus. This points to the 
possibility of unifying the results and generalizing to arbitrary complex functions. 
An assumption of analyticity is required; technically, this is stronger than smooth- 
ness but in practice one is never satisfied without the other. Assuming analyticity, 
derivations in the Fourier and Laplace cases may indeed by unified via a hybrid 
approach. In one variable, this is carried out all the time under various names such 
as "steepest descent", "stationary phase" or "saddle point". Writing 1(A) as a 
complex contour integral, the critical point will be a saddle point in for the real 
part of the phase function; the contour may then be re-oriented to pass through the 
saddle in the direction of steepest descent of —(j), converting the integral into one 
of Fourier type and explaining why the series are nearly identical. This is carried 
out, for example, in [dB81] or [BH86, Chapter 7]. 

In more than one variable, complex phase results generalizing (1.3) and (1.5) 
are strangely absent from the literature. They are not treated in [Ste93; Hen91; 
BH86; Won89; Bre94; dB81]. Wo have found only two significant treatments of 
complex phases. One is Section 7.7 of [H6r83]. Here, notions from several complex 
variable theory are used to validate integration by parts, leading to asymptotic 
expansions in the case when A is smooth and compactly supported. This work has 
been used by [RW08] to produce explicit higher order terms in various expansions. 
As stated, Hormander's method does not allow integration over manifolds with 
boundary. The second treatment is [Fed89], which uses similar methodology to 
ours (complex contour deformation). The cited work omits the details, referring to 
the original article [Fed77] in Russian. Regions of integration with boundaries are 
allowed here, but only if Rc{(/)} is strictly maximized on the interior (our preliminary 
case, discussed in Section 4). In applications, work requiring high-dimensional 
saddle integration with complex phases is usually done by hand; see, e.g., [Wor04] 
or [BFSSOO]. 

There are several possible reasons for this gap. One is that fewer people are 
aware of the corresponding techniques in several complex variables. Results such 
as Cauchy's integral theorem in several variables and Stokes' theorem for holomor- 
phic d- forms in C^, while known to first-year graduate students in several complex 
variables, are not always known to those applying these techniques. The few math- 
ematicians who do use these techniques to evaluate integrals [BH91; Lic91] are not 
unusually interested in anything as pedestrian as the generalization of standard 
saddle point integral evaluation. Secondly, even if one knows the techniques, mov- 
ing contours in higher dimensions is a tricky business. The geometry can be much 
more complicated. When the real part of the phase vanishes at points that are 
not critical, further care must be taken in how the contour is moved. Finally, the 
need for such a generalization, illuminating though it may be, may have been too 
infrequent. 

The organization of rest of the paper is as follows. In the next section we 
give some definitions having to do with stratified spaces. We then state our main 
result. Theorem 2.3. Section 3 records some easy computations in the case (p is the 
standard phase function S'(x) := J2j=i Section 4 handles the general case under 
the assumption that the real part of (p has a strict minimum at the critical point. 
Theorem 2.3 is proved in Section 5. Section 6 gives an application from [PW04] to 
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the estimation of coefficients of a bivariate generating function. Finally, Section 7 
discusses further research directions. 

2. Notation and statement of results 

Stratified spaces. Because of the useful properties ensuing from the defini- 
tion, wc will use Whitnc'v wtratificd spaces as our chains of integration. Aside from 
these useful properties, the details of the definition need not concern us, though 
for completeness we give a precise definition. Let / be a finite partially ordered set 
and define an /-decomposition of a topological space Z to be a partition of Z into 
a disjoint union of sets {Sq, : a G /} such that 

S'„n5^7^0 ^ So,^~S^i ^ a<l3. 

Definition 2.1 (Whitney stratification). Let Z be a closed subset of a smooth 
manifold A4. A Whitney stratification of Z is an I -decomposition such that 
(i) Each Sa is a manifold in R". 

(a) If a < (3, if the sequences {xt G S/)} and {y, e Sa} both converge to y £ Sa, if 
the lines U = xiyl converge to a line I and the tangent planes Txi{Sp) converge 
to a plane T of some dimension, then both I and Ty{Sa) are contained in T. 

For example, any manifold is a Whitney stratified space with one stratum; any 
manifold with boundary is a Whitney stratified space with two strata, one being 
the interior and one the boundary; a fc-simplex is a Whitney stratified space whose 
strata are all its faces. Whitney stratified spaces are closed under products and the 
set of products of strata will stratify the product. Every algebraic variety admits a 
Whitney stratification, although the singular locus filtration may be too coarse to 
be a Whitney stratification. 

Critical points. Associated with the definition of a stratification is the strat- 
ified notion of a critical point. Observe that under this definition, any zero- 
dimensional stratum of yVI is a critical point of A4. 

Definition 2.2 (smooth functions and their critical points). Say that a func- 
tion (j) : M. ^ €. on a stratified space A4 is smooth if it is smooth when restricted 

to each stratum. A point p € M is said to be critical for the smooth function <j) if 
and only if the restriction d<j}\s vanishes, where S is the straium containing p. 

Let C be a real analytic, d-dimensional stratified space. This means that 
each stratum 5 is a subset of and each of the chart maps ijj from a neighborhood 
of the origin in R'^ to some fc-dimcnsional stratum 5 C C'' is analytic (the coordinate 
functions are convergent power series) with a nonsingular differential. It follows 
that V may be extended to a holomorphic map on a neighborhood of the origin 
in C'^, whose range we denote by 5 (g) C. Choosing a small enough neighborhood, 
we may arrange for 5 ® C to be a complex fc- manifold embedded in C*. We say 
a function / : — > C is analytic if it has a convergent power series expansion in 
a neighborhood of every point; an analytic function on M may be extended to a 
complex analytic function on a neighborhood of in (g) C := Uae/ ® 
Because we are interested in the integrals of rf- forms over M., there is no loss of 
generality in assuming that A4 is contained in the closure of its rf-dimcnsional 
strata, whence AI (g) C is a neighborhood of M. in C^. Real d-manifolds in 
are not naturally oriented, so we must assume that an orientation is given for 
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each (i-stratuni of A4, meaning that the chart maps from R'' to M must preserve 
orientations. The critical point p € is said to be quadratically nondegenerate 
if p is in a d-dimensional stratum and Hip) is nonsingular, where 



nip) := 



dxi dxj / 



is the Hessian matrix for the function tj) of Definition 2.2 on a neighborhood of p in 

Finally, we define a point of M. to be stationary if it is critical and Bje{(f)} 
vanishes there. 

Results. The main result of the paper is an asymptotic expansion for X(A). 

Theorem 2.3. Let M. he a com,pact, real analytic stratified space of dimension 
d embedded in C^. Let A and tp be analytic functions on a neighborhood of M. and 
suppose Re{(l)} > on M.. Let G be the subset of stationary points of (p on M and 
assume that G is finite and each stationary point is quadratically nondegenerate. 
Then the integral 2{X) := Jj^e~^'^^'^^A{x) dx has an asymptotic expansion 

oo 



// A is nonzero at some point of G then the leading term is given by 
(2.1) Co = {2Tr)-'^/^^A{x)uj{x,X){detH{x)y^^^ . 

where uj{x, A) is the unit complex number e~^'^^^h The choice of sign of the square 
root on the right-hand side of (2.1) is determined by choosing any analytic chart 
map ijj for a neighborhood of x and defining 

(detH(a;))-i/2 ._ (detH((/) o V))"^/^ det 

where is the Jacobian matrix for ip at x and the —1/2 power of the determinant 
on the right is the product of the inverses of principal square roots of the eigenvalues. 

Corollary 2.4. Assume the hypotheses of Theorem 2.3 except that some of 
the points of G may be in strata of dimension d — 1 with a neighborhood in M. 

diffeom,orphic to a half space in . Then the same conclusion holds except that the 
summand in (2.1) corresponding to such a point x must be multiplied by 1/2. 

3. Preliminary results for the standard phase function 

For X e W^, let S'(x) := denote the standard quadratic form. We 

begin with a couple of results on integrals of Laplace type where the phase function 
is the standard quadratic and A is a monomial. In one dimension, 



Proposition 3.1. 



/. 



OO I 

e-^-' dx = ^^—^ ;,-l/2-n/2 

(n/2)!2"\/2^ 



if n is even, while the integral is zero if n is odd. 
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Proof. For odd n the result is obvious from the fact that the integrand is 
an odd function. For n = the result is just the standard Gaussian integral. By 
induction, assume now the result for n — 2. Integrating by parts to get the second 
line, we have 

„n~l 



n — 1 
2A 

n-1 1 (n-2)! 



x"' e dx 

^l/2-n/2) 



2A ^(n/2-l)!2"-i 

by the induction hypothesis. This is equal to A-V2-"/2(27r)-i/2n! / ((n/2)! 2"), 
completing the induction. □ 

Corollary 3.2 (monomial integral). Let r be any d-vector of nonnegative 
integers and let x'' denote the monomial x^^ • • • x^f . Then 



I 

Js. 



where 



(3.1) /3r := (27r)-'^/2 -Q 



if all the components rj are even, and zero otherwise. 
Proof. The integral factors into 

d r 

n 



■oo 



□ 



To integrate term by term over a Taylor series for A, we need the following 
estimate. 

Lemma 3.3 (big-0 estimate). Let A he any smooth function satisfying A{x) = 
Odxl*") at the origin. Then the integral of A{x)e~^^'^^^ over any compact set K 
may he hounded from ahove as 



I A(x)e-^^(^) dx = 0(A 



-(d+r)/2 



Proof. Because K is compact and ^(x) = 0(|x|'') at the origin, it follows 
that there is some constant c for which |^(x)| < c|x|'' on all oi K. Let Kq denote 
the intersection of K with the ball |x| < A~^/^ and for n > 1 let Kn denote the 
intersection of K with the shell 2"-^-^/^ < |x| < 2"A-V2. On Kq we have 

|A(x)| < cA-''/^ 

while trivially 

/ e-^^Wdx< / dx<CdA-'^/2. 

JKn JKn 
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Thus 



/ ^(x)e-^^(^) dx 

Jko 



For n>l, on K C\ Kn, we have the upper bounds 



-AS(x) 



< 



I 



dx 



Letting c" 



< 2™ Cd X-^/'^ . 
< oo, we may sum to find that 



n=0 

proving the lemma. 



oc ^ 

V / A(x)e-^^(^) rfx 



< (c'+c")A-('^+<^)/^ 



□ 



It is now easy to compute a series for 1(A) in the case where is the standard 
quadratic and the integral is over a neighborhood of the origin in M"^. 

Theorem 3.4 (standard phase). Let A{x) = Orx'' be any analytic function 
defined on a neighborhood J\f of the origin in M**. Let 



(3.2) 

Then 



J(A) := / A(x)e-^^W(ix. 



i{x) ~ 



^-(n+<i)/2 



_|r|=n 

as an asymptotic series expansion in decreasing powers of X, with (5r as in (3.1). 
Proof. Write A{x) as a power series up to degree A'' plus a remainder term: 

(E E "rxM +i?(x) 

where i?(x) = 0(|x|^+^). Using Corollary 3.2 to integrate all the monomial terms 
and Lemma 3.3 to bound the integral of R{x)e~^^^^^ shows that 

N 

1(A) = ^ ^ ar/?rA-('^+")/2 + o(A-('^+")/2-i) 

n=0 |r|=n 

which proves the asymptotic expansion. 



□ 



4. The case of a strict minimum 



In this section, we continue to integrate over a neighborhood of the origin in 
M**, but we generalize to any analytic phase function <p with the restriction that 
the real part of 4> have a strict minimum at the origin. The assumption of a strict 
minimum localizes the integral to the origin, so the only tricky aspects are keeping 
track of the sign (Lemma 4.3) and being rigorous about moving the contour. 



8 



ROBIN PEMANTLE AND MARK C. WILSON 



Theorem 4.1. Let A and (p be complex-valued analytic functions on a compact 
neighborhood Af of the origin in M'' and suppose that the real part of cj) is nonneg- 
ative, vanishing only at the origin. Suppose that the Hessian matrix Ti. of (j) at the 
origin is nonsingular. Denoting 1(A) := A{-x)e~'^'^^^^ dx, there is an asymptotic 
expansion 

i>0 

where 

and the choice of sign is defined by taking the product of the principal square roots 
of the eigenvalues ofH. 

The proof is essentially a reduction to the case of standard phase. The key 
is the well known Morse Lemma. The proof given in [Ste93, VIII:2.3.2] is for the 
smooth category and purely real or imaginary phase but extends without significant 
change to complex values and the analytic category. For completeness, we include 
the adapted proof. 

Lemma 4.2 (Morse lemma). There is a bi-holomorphic change of variables 
X = V(y) such that (f){tp{y)) = S{y) := Y^^^i y]- The differential = dip{0) will 
satisfy (det J^)^ = (det ^H)-^. 

Proof. Addressing the second conclusion first, we recall how the Hessian ma- 
trix behaves under a change of variables. If ip : ^ is bi-holomorphic on a 
neighborhood of x and if (p has vanishing gradient at and Hessian matrix H 
there, then the Hessian matrix ?i of (/) o at x is given by 

where is the Jacobian matrix for tp at x. The standard form S has Hessian 
matrix equal to twice the identity, hence any function ip satisfying (po'ip = S must 
satisfy 

2 Id = Ti, J,p . 

Dividing by two and taking determinants yields |J^pdet(i?i) = 1, proving the 
second conclusion. 

To prove the change of variables, the first step is to write 

d 

where the functions (pj^k = (pkj are analytic and satisfy (j)j^k{0) = (l/2)Hj,fe- It 
is obvious from a formal power series viewpoint that this can be done because the 
summand XjXkCpj^^ can be any power series with coefficients indexed by the orthant 
{r : r > (5j -|- 5^}; these orthants cover {r : |r| > 2}, so we may obtain any function 
^ vanishing to order two; matching coefficients on the terms of order precisely two 
shows that (^j,fc(0) = (1/2) Wj,^- 

More constructively, we may give a formula for (j)j^k. There is plenty of freedom, 
but a convenient choice is to let Or denote the coefficient of x'' in 0(x) and to take 

XkXk(pj,k{'^) ■= } ^ , ,n I T\- OlrX . 

^-^ r ( r — 1) 

|r|>2 I I 
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For fixed r, it is easy to check that 



l<j,k<d 



M(|r|-1) 



whence (t> = Yl XjXk4>j,k- Alternatively, the following analytic computation from [Ste93] 

verifies that (p = kXjXk4>j,k- Any function / vanishing at zero satisfies f{t) = 

— s)f'{s) ds, as may be seen by integrating by parts (take g{s) = —(1 — s)). 
Fix X and apply this with f{t) = {d/dt)(j){tx) to obtain 

"/•(x) = ^ ^(/-(tx) dt = t)^Ht^) dt . 



The multivariate chain rule gives 



— (/)(fx) =}_^XjXk 



dxjdxk 



(tx); 



plug in = ttrx"" and integrate term by term using — t)t'^ dt = 

to see that cj) = Y^j.k ^3^'k4>3.k- 

The second step is an induction. Suppose first that '/>jj (0) ^ for all j. The 
function 0^ ^ and a branch of the function (^^'-^ are analytic in a neighborhood of 
the origin. Set 

,1/2 , yk4>i,k 

fe>i ^ . 

Expanding, we find that the terms of of total degree at most one in the terms 
X2, - ■ ■ ,Xci match those of (f) and therefore, 

(4.1) (p{x) = + ^ XjXkhj^k 

j,k>2 

for some analytic functions hj^k satisfying hj^k{0) = {1/2) Hj^k- Similarly, if 

r-l 



(x) = J2yj+ ^3^khj,k 



then setting 



gives 



yk(f>r,k 



k>r 



j=l j,k>r+l 

for some analytic functions hj^k still satisfying /ij,fe(0) = {l/2)Hj^k- By induction, 
we arrive at <^(x) = J2'j=i v'ji finishing the proof of the Morse Lemma in the case 



where each Ti.j^j is nonzero. 



Finally, if some Hj^k = 0, because Ti. is nonsingular we may always find some 

unitary map U such that the Hessian U'^ HU oi (j> o U has no vanishing diagonal 
entries. We know there is a tpo such that {(j)oU) otpo = S, and taking tp = U o tpQ 
finishes the proof in this case. □ 
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Proof of Theorem 4.1: The power series allows us to extend </) to a neighborhood 
of the origin in C*. Under the change of variables ^jJ from the previous lemma, we 
see that 

J(A) = / ^oV'(y)e-^^('')(det#(y))rfy 
:= / i(y)e-^^('^) dy 

for some analytic function A, where C is a neighborhood of the origin in M". We 
need to check that wc can move the chain tp^^C of integration back to the real 
plane. 

Let h{z) Rc{S'(z)}. The chain C := V^H^) lies in the region {z e : 
/i(z) > 0} except when z ~ 0, and in particular, h > e > on dC . Let 

H{z, t) := Re{z} + {l-t)i Im{z} . 

In other words, i? is a homotopy from the identity map to the map tt projecting 

out the imaginary part of the vector z. For any chain cr, the homotopy H induces 
a chain homotopy, H{a) supported on the image of the support of a under the 
homotopy H and satisfying 

dH{a) = CT - TTCT + H{da) . 

With <j = C, observing that S{H{z,t)) > S{z), we see there is a l)-chain V 
with 

dV^C- nC + C" 

and C" supported on {h> e}. Stokes Theorem tells us that for any holomorphic 
d-form w, 



/ " = / 

JdV Jv 



(Lj = 



and consequently, that 



iJ = U) + CO . 

Jc JwC JC" 

When Lo = Ae~^^ dy, the integral over C" is 0{e~^^), giving 
J(A) = / A{y)e-^^^y^ dy + 0(6"^^) . 

JttC 

Up to sign, the chain ttC" is a disk in R'' with the standard orientation plus 
something supported in {h> e}. To see this, note that tt maps any real d-manifold 
in diffeomorphically to M.'^ wherever the tangent space is transverse to the imag- 
inary subspace. The tangent space to the support of C at the origin is transverse 
to the imaginary subspace because S > on C, whereas the imaginary subspace 
is precisely the negative c?-spacc of the indcx-c? form S. The tangent space varies 
continuously, so in a neighborhood of the origin, tt is a diffeomorphism. Observing 
that i(0) = A(0) det(rfV'(0)) = ^(0)(det W)-^^ and using Theorem 3.4 finishes 
the proof up to the choice of sign of the square root. 

The map dn o dip~^{0) maps the standard basis of M'' to another basis for R'^. 
Verifying the sign choice is equivalent to showing that this second basis is positively 
oriented if and only if dct(dip(0)) is the product of the principal square roots of the 
eigenvalues of H (it must be either this or its negative). Thus we will be finished 
by applying the following lemma (with a = V'~^)- 
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Lemma 4.3. Let C fee the set {z : Rc{S'(z)} > 0}. Pick any a G GLd{C) 
mapping K*^ into W and let M :— a he the matrix representing S o a. Let 
TT : ^ M."^ be projection onto the real part. Then iroa is orientation preserving on 
M'^ if and only if det a is the product of the principal square roots of the eigenvalues 
of M (rather than the negative of this). 

Proof. First suppose a e GLa{R). Then M has positive eigenvalues, so the 
product of their principal square roots is positive. The map tt is the identity on 
M'' so the statement boils down to saying that a preserves orientation if and only 
if it has positive determinant, which is true by definition. In the general case, let 
at := TTt o a, where 7rt(z) = Re{z} + (1 — t)lm{z}. As we saw in the previous 
proof, 7rt(IR.'') C W for all < t < 1, whence Mt := afat has eigenvalues with 
nonnegative real parts. The product of the principal square roots of the eigenvalues 
is a continuous function on the set of nonsingular matrices with no negative real 
eigenvalues. The determinant of at is a continuous function of t, and we have seen 
it agrees with the product of principal square roots of eigenvalues of Mt when t = 1 
(the real case), so by continuity, this is the correct sign choice for all < t < 1: 
taking t = proves the lemma. □ 

For later use, we record one easy corollary of Theorem 4.1. 

Corollary 4.4. Assume the hypotheses of Theorem 4-1 and let M' he the 
intersection of M with a region diffeomorphic to a half space through the origin. Lf 
A{Q) ^ then 

:= / ^(x)e-^^W dxr.^ X-"'^ 
where cq is the same as in the conclusion of Theorem 4-1- 

Proof. Under the change of variables and the projection tt, this region 
maps to a region J\f" diffeomorphic to a halfspace with the origin on the boundary. 
Changing variables by y = A~^/^x and writing M\ for X^/^M", we have 

T'{X)=X-^'^ I Ax{y)e-^^^Uy 

where ^^(y) = {Aotp){X^^/^y). The function Ax converges to A{0) pointwise but 
also in where fi is the Gaussian measure e"'^^^) dx. Also, the regions Af\ 

converge to a halfspace H in the sense that their indicators l^v^ converge to 1h in 
L^(^). Thus A\lj^^ converges to A{0)1h in and unravelling this statement 

we see that 

/ Ax{y)e-'^y^ dy ^ [ A(0)e-^(y)rfy. 

The last quantity is equal to Co/2, showing that X'^/'^I'{X) — > Co/2 and finishing the 
proof. □ 

5. Proofs of main results 

Theorem 2.3 differs from Theorem 4.1 in several ways. The most important 
is that the set where ^(j) vanishes may extend to the boundary of the region of 
integration. This precludes the use of the easy deformation tt because C" is no 
longer supported on {h > e}. Consequently, some work is required to construct a 
suitable deformation. We do so via notions from stratified Morse theory [GM88] . 
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Tangent vector fields. If x is a point of the stratum S of the stratified space 
M, let Tx(A1) denote the tangent space to S at x. Because M is embedded in 
C^, the tangent spaces may all be identified as subspaces of C''. Thus we have a 
notion of the tangent bundle TM, a section of which is simply a vector field / on 
M C C''- such that /(x) e Tx(A1) for all x. A consequence of the two Whitney 
conditions is the local product structure of a stratified space: a point p in a A;- 
dimensional stratum S of a, stratified space A4 has a neighborhood in which A4 is 
homeomorphic to some product S xX. According to [GM88], a proof may be found 
in mimoegraphed notes of Mather from 1970; it is based on Thom's Isotopy Lemma 
which takes up fifty pages of the same mimeographed notes. The next lemma is 
the only place where we use this (or any) consequence of Whitney stratification. 

Lemma 5.1. Let f be a smooth section of the tangent bundle to S, that is 
f{s) £ Ts{S) for s e S. Then each s E S has a neighborhood in A4 on which f 
may be extended to a smooth section of the tangent bundle. 

Proof. Parametrize M locally by 5 x A and extend / by f{s, x) := f{s). □ 

Lemma 5.2 (vector field near a non-critical point). Let x be a point of the 
stratum S of the stratified space M. and suppose x is not critical for the function (j). 
Then there is a vector v G T:x.{S® M.) such that Re{(i(/)(v)} > x. Furthermore, 
there is a continuous section f of the tangent bundle in a neighborhood M o/x such 
that Re{rf(?!)(/(y))} >Q at every y GJ\f. 

Proof. By non-criticality of x, there is a w e Tx(5) with d(p{w) = u 7^ at 
X. Multiply w componentwise by u to obtain v with Rc{(i0(v)} > at x. Use any 
chart map for 5 (g) C near x to give a locally trivial coordinatization for the tangent 
bundle and define a section / to be the constant vector v; then Re{c?<?!)(/(y))} > on 
some sufficiently small neighborhood of x in S. Finally, extend to a neighborhood 
of X in by Lemma 5.1. □ 

Although we are working in the analytic category, the chains of integration are 
topological objects, for which we may use C°° methods (in what follows, even 
methods will do). In particular, a partition of unity argument enhances the local 
result above to a global result. 

Lemma 5.3 (global vector field, in the absence of critical points). Let A4 be 
a compact stratified space and (j) a smooth function on A4 with no critical points. 
Then there is a global section f of the tangent bundle of M. such that the real part 
of d(j){f) is everywhere positive. 

Proof. For each point x e Ad, let /x be a section as in the conclusion of 
Lemma 5.2, on a neighborhood ?7x- Cover the compact space M by finitely many 
sets {C/x : X e F} and let {V'x : x e .F} be a smooth partition of unity subordinate 
to this finite cover. Define 

/(y) = E^-(y)/-(y)- 

Then / is smooth; it is a section of the tangent bundle because each tangent space 

is linearly closed; the real part of d(f){f{y)) is positive because we took a convex 
combination in which each contribution was nonnegative and at least one was pos- 
itive. □ 
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Another partition argument gives our final version of this result. 

Lemma 5.4 (global vector field, vanishing only at critical points). Let A4 be a 

compact stratified space and a smooth function on Ai with finitely many critical 
points. Then there is a global section f of the tangent bundle of M. such that the 
real part of d(j){f) is nonnegative and vanishes only when y is a critical point. 

Proof. Let A4e be the compact stratified space resulting in the removal of an 
e-ball around each critical point of cj). Let be a vector field as in the conclusion of 
Lemma 5.4 with A4 replaced by A4e. Let c„ be a positive real number, small enough 
so that the magnitudes of all partial derivatives of c„/i/„ of order up to n are at 
most 2'". In the topology of uniform convergence of derivatives of bounded order, 
the series '}2n Cnfn converges to a vector field / with the required properties. □ 

Proof of Theorem 2.3. Let / be a tangent vector field along which (j) in- 
creases away from critical points, as given by Lemma 5.4. Such a field gives rise to 
a differential flow, which, informally, is the solution to dp/dt = /(p). To be more 
formal, let x be a point in a stratum S of M. Via a chart map in a neighborhood 
of X, we solve the ODE d^{t)/dt = with initial condition $(0) = x, ob- 

taining a trajectory $ on some interval [0, Cx] that is supported on S. Doing this 
simultaneously for all x e results in a map 



with <i?(x, t) remaining in 5 (g) C when x is in the stratum S. The map $ satisfies 
$(x,0) = X and {d/dt)^{x,t) = /($(x,t)). The fact that this may be defined 
up to time e for some e > is a consequence of the fact that the vector field / 
is bounded and that a small neighborhood of M in (8) C is embedded in C^. 
Because / is smooth and bounded, for sufficiently small e the map x 1— > $(x, e) is 
a diffeomorphism. 

The flow reduces the real part of 4> everywhere except the critical points which 
are rest points. Consequently, it deflncs a liomotopy H{K,t) := $(x, et) between C 
and a chain C on which the minima of the real part of <f> occur precisely on the set G. 
Recall that H induces a chain liomotopy Ch with dCn = C — C + dC x where a is 
a standard 1-simplex. Let lu denote the lioloiiiorphic d-form A(z) exp(— A0(z)) dz. 
Because a; is a holomorphic d-form in C^, we have du) = 0. Now, by Stokes' 
Theorem, 



The chain dC x a is supported on a finite union of spaces S^C where S* is a stratum 
of dimension at most d — 1. The integral of a holomorphic d-form vanishes over 
such a chain. Therefore, the last term on the right drops out and we have 



Outside of a neighborhood of G the magnitude of the integrand is exponentially 
small, so we have shown that there are rf-chains Cx supported on arbitrarily small 



$ : X [0, e] ^ C 



id 
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neighborhoods A/'(x) of each x e G such that 



(5.1) J(A) - E / 

is exponentially small. To finish that proof, we need only show that each ui has 
an asymptotic series in decreasing powers of A whose leading term, when A{x) ^ 0, 

is given by 

(5.2) co(x) = (27r)-'^/2^(x)e^*('')(detH(x))-i/2 

The d-chain Cx may by parametrized by a map V'x : -B -^(x)) mapping the origin 
to X, where B is the open unit ball in M''. By the chain rule, 

w = / [Ao exp{—X[(j) o ^{x}]) detd-)p{'x.) dii. 

JB 

The real part of the analytic phase function <p o ijj has a strict minimum at the 
origin, so we may apply Theorem 4.1. We obtain an asymptotic expansion whose 
first term is 

(5.3) (27rA)-'*/2[yl o V](0) det rfV(0)(det M^)"^/^ 

where Mx is the Hessian matrix of the function (poip. The term [Ao^/;] (0) is equal to 
A{x). The Hessian matrix of at the origin is given by Mx = d'ip{0) H(x) dtp{0). 
Thus 

detMx = (det#(0))^detW(x) 
and plugging into (5.3) yields (5.2). □ 



Proof of Corollary 2.4: Lemma 5.4 does not require the critical points to be in 
the interior, so the argument leading up to (5.1) is still valid. For those points 
X in a (d — l)-dimensional stratum, use Corollary 2.4 in place of Theorem 4.1 to 
obtain (5.2) with an extra factor of 1/2. □ 

Remark. The reason we do not continue with a litany of special geometries 
(quarter-spaces, octants, and so forth) is that the case of a half space is somewhat 
special. The differential of the change of variables at the origin is a nonsingular 
map, which must send half-spaces to half-spaces, though it will in general alter 
angles of any smaller cone. 



6. Examples 

The simplest multidimensional application of our results is a computation 
from [PW04]. The purpose is to estimate coefficients of a class of bivariate gener- 
ating functions whose denominator is the product of two smooth divisors. We give 
only a brief summary of how one arrives at (6.1) from a problem involving gener- 
ating functions; a complete explanation of this can be found in [PW04, Section 4]. 
Note, however, that the mathematics of the integral is not contained in that paper, 
which instead refers to an earlier draft of this one! 

Let v\,V2 be distinct analytic functions of z with v\{l) =^2(1) = l,Q ^ v'i{l) ^ 
^2(1) 7^ 9,nd such that each |i-'i(z)| attains its maximum on \z\ = 1 only at z = 1. 
For example, the last condition is satisfied by any pair of aperiodic power series 
with nonnegative coeSicients and radius of convergence greater than 1. 
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Consider the generating function F{z,w) — 1/H{z,w), where H{z,w) = (1 — 
wvi{z)){l — wv2{z)). The two branches of the curve H = intersect only at 
(1, 1) and this intersection is transverse. The Maclaurin coefficients of F{z,iu) = 
J2r s O'TsZ^W' are given by the Cauchy integral formula 

1 dwdz 
"'''^ " (27ri)2 - wvi{z)){l - wv2{z)) 

where the integral is taken over a product of circles centred at (0, 0) and of suffi- 
ciently small radii. 

Pushing the contour out to \z\ = = 1 — ewe obtain the same formula, 

since F is still analytic inside the product of disks boimdcd by these latter circles. 
Pushing the ui-contour out to |w| = 1 + e, using the residue formula on the inner 
integral and observing that the integral over \w\ = 1 + e is exponentially decaying 
as s ^ oo, we see that 

-Rs{z) 



■ dz 



where ~ means that the difference is exponentially decaying as s ^ oo and Rs{z) 
denotes the sum of residues of w i— > w~^^~^^^ F{z,w) at the roots w = l/vi{z), 

{1,2}. 

The residue sum Eg (z) can be rewritten in terms of an integral via 



-Rs{z) = (s + 1) / [(1 - p)vi + PV2Y dp 
Jo 

[ ^-('■+1) [' [{l-p)vi{z) +pv2{z)Y dpdz. 

J\z\ = l Jo 



and so we have 

In order to cast this into our standard framework, we need to be able to define a 
branch of the logarithm of (1 — p)vi{z) + pv2{z). We do this by localizing on the 
circle \z\ = 1 to a sufficiently small neighbourhood of the point z = 1. This is 
possible since the integrand decays exponentially away from z = 1, by hypotheses 
on the Vi, and we shall show that the integral near z = 1 decays only polynomially. 
The substitution z = e'^ converts this to an integral 

(6.1) a„«^-^/ / e-'"l'^P^'^A{p,z)dpdt 

Sti" Jm Jo 

where (f>{p, t) = ir0/s + \og [(1 — p)vi{e^*) + 15^2 (r*')] , A{p, t) = 1, and A/" is a closed 
interval centred at 0. To compute asymptotics in the direction r/s = k, for fixed 
K > 0, we can consider (f) to be independent of r and s. 

We now asymptotically evaluate (6.1) using Theorem 2.3. We can rewrite the 
iterated integral as a single integral over the stratified space M = Af x [0, 1]. The 
phase (f) has nonnegative real part and this fits into our framework. There is a 
single stationary point, at (p, z) = (1/2, 0) (note that Rc{(/)} is zero for all (p, 0), so 
Theorem 4.1 does not suffice). This critical point is quadratically nondegenerate 
and direct computation using Theorem 2.3 yields 

(6-2) ars = I ^ +0{s-') 

bi(l) -^'2(1)1 

as s ^ 00 with K fixed. By keeping track of error terms more explicitly, it is easily 
shown that this approximation is uniform in k provided k stays in a compact subset 
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of the open interval formed by w^(l),f2(l) (it foUows from our assumptions that 
these numbers are positive real — see [PW04] for more details). This means that 
ttrs is asymptotically constant in any compact subcone of directions away from the 
boundary formed by the Hues Ki = v'^{l). 

This example, and in fact a number of cases in [PW04], can also be solved 
using iterated residues. This is carried out in [BP04]. Iterated residues have the 
advantage of showing that the 0{s~^) term decays exponentially, but the disad- 
vantage that they do not give any results when k approaches the boundary. The 
present methods do give boundary results. Corollary 2.4 shows that ttrs converges 
to one half the right-hand side of (6.2) when (r, s) — > oo with r/s = K\ + 0(1), and 
a small extension obtains a Gaussian limit: letting $ denote the standard normal 
cumulative distribution function, we have 

" - +0(a-) 



K(i)-^'^(i)l 

when r, s ^ oo with (r/s — Ki)/s^/'^ u. 



7. Further comments 

We have not emphasized explicit formulae for the higher order terms, giving 
an equation such as (2.1) only for the leading term in the case where A{Qi) ^ 0. 
However, our results establish the validity of existing computations of higher order 
terms under our more general hypotheses. 

To elaborate, we prove Theorem 2.3 by first constructing a change of variables 
X I— > $(a;, e) homotopic to the identity under which the minimum of Re{(/>} at 
is strict, and then changing variables, again homotopically to the identity, to the 
standard form. The composition ij} of these two maps is homotopic to the identity 
but is far from explicitly given: while the second map is constructed by an explic- 
itly defined Morse function, the; first deformation is the solution to a differential 
equation and is not particularly explicit. 

In [Hor90] , Hormander derives such an explicit formula (assuming smoothness) 
for integrals of our type where M —W'' and A has compact support. The formula is 
indeed rewritten and used in [RW08] to compute higher order terms for generating 
function applications, in which more restrictive hypotheses preclude the vanishing 
of Re{0} on a curve reaching the boundary of the chain of integration. Their 
methods, while not covering the cases of interest here, do have the virtue of dealing 
with the change of variables tj) only through the equation S = (j)0'tp. In particular, 
the derivatives of ip arising in the computation of the new amplitude function 
(Aotp) dctdip can be computed by implicitly differentiating the equation S = (fx^^jj. 
Having found at least one such ip homotopic to the identity, we are now free to 
replicate the computations of [RW08] under our more general hypotheses, as follows. 

In the case of standard phase, the coefficient of A"*^""'"'')/^ is given by (provided 
all Ti are even) 

ar/3r 

|r|=n 

where Or is the Maclaurin coefficient of A corresponding to the monomial r and 
(3r is the constant defined in Corollary 3.2. Note that n must be even for this 
coefficient to be nonzero, so we write n = 2k. The differential operator . . . d^^ 
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when applied to A and evaluated at yields precisely Yii Tilar- Thus the operator 

a^'^i . . . ^ 



E 

|r|=fc 



. . . r^! 



applied to A and evaluated at zero yields the coefficient wc seek. 

After the Morse lemma is applied using the change of variables S = o i/i, 
we need to apply the displayed operator to the new amplitude (A o tp) det d-ip. The 
resulting expression evaluated at x can be computed directly via the rules of Leibniz 
and Faa di Bruno. Evaluating at x simplifies some terms, and, as mentioned above, 
derivatives of (A o ip) det dip may be computed without explicitly specifying ^p. 

As a relatively simple example, consider the case k = 1 and d = 1. The 
differential operator reduces to jd"^ where d denotes differentiation with respect to 
the variable x. Applying this to {A o ^p) det dip yields (with superscripts denoting 
the order of derivatives and arguments suppressed) 



^(2) 



The defining equation S — cpotp can be differentiated to yield the system 

2 = </,(2) L(l)l' + ,^(1)^(2) 



= </)(3) 



= 0(4) ^(1) +60(3) V-^i) V'^2)+40(2)V'(i)t/'(3)+30(2) V(') +0(1)^'^^). 

Evaluating these at the point in question, we see that the terms with highest 
derivatives of \p vanish in each equation. The system is triangular and can be solved 
explicitly to obtain 




v>(3) 



3 [0(2)]' 

[50(-^)]' -30(2)0(4) 
3V2 [0(2)] 



7/2 



Putting these together with the expression for the derivative of {Aoxp) det dtp above 
yields an expression for the A~^/2 term in the integral that is a rational function 

with denominator [0(2)] '^^^ and numerator a polynomial in the derivatives of A up 
to order 2 and to order 4. In summary, the results of this paper show that the 
computational apparatus and formulae for higher order terms given in [RW08] hold 
in the case of complex phase functions integrated over stratified spaces. 
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